Self-Organization of Vortex Length Distribution in Quantum Turbulence: 
An Approach from the Barabasi- Albert Model 
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The energy spectrum of quantum turbulence obeys Kolmogorov's law. The vortex length distri- 
bution (VLD), meaning the size distribution of the vortices, in Kolmogorov quantum turbulence 
also obeys a power law. We propose here an innovative idea to study the origin of the power law of 
the VLD. The nature of quantized vortices allows one to describe the decay of quantum turbulence 
with a simple model that is similar to the Barabasi- Albert model of large networks. We show here 
that such a model can reproduce the power law of the VLD fairly well. 

PACS numbers: 03.75.Lm, 07.05.Mh, 05.40.-a, 67.40.Vs, 47.37.+q, 89.20.Hh, 84.35.+i 
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Scale-free (scale invariant) systems are ubiquitous. Ex- 
amples of scale-free systems include the World Wide 
Web (WWW), whose vertices are HTML documents con- 
nected by hyperlinks pointing from one page to another, 
and social networks such as a movie collaboration net- 
work of movie actors in which each actor is represented 
by a vertex, two actors being connected if they were cast 
together in a same movie. Regardless of the system and 
the identity of its constituents, the probability P(k) that 
a vertex in the network connects with other k vertices 
decays as a power law, following P(k) oc fc -7 Q. 

Fluid turbulence has a similar scaling property. Theo- 
retical and experimental studies have shown that classical 
turbulence has the scale invariance in its energy spectrum 
E{k) = 47rfc 2 |w(fc)| 2 , where v(k) is the Fourier transform 
of fluid velocity v(r) 0. This energy spectrum obeys the 
Kolmogorov power law E{k) oc /c~ 5 / 3 , which means that 
classical turbulence has scaling property in wave num- 
ber space. Some numerical simulations recently revealed 
that decaying quantum turbulence at zero temperature, 
in which quantized vortex lines of various length form an 
irregular tangle Q , is also a scale invariant system 0, |j| ; 
quantized vortices are topological defects of order param- 
eter in a superfluid. These authors calculated the energy 
spectrum of the velocity field and found that it obeys 
the Kolmogorov power law E(k) oc k~ 5 ^ 3 as in classical 
turbulence. The Kolmogorov law is a scaling property in 
wave number space, but should be closely related with 
the self-similarity of the turbulent velocity field in real 
space. So we ask the following question: how can the 
scale invariance in wave number space be translated into 
the scale invariance in real space? A possible answer was 
proposed by Araki et al. 4] . They calculated the vortex 
length distribution (VLD) n(l), where n(l)Al represents 
the number of vortices in Kolmogorov quantum turbu- 
lence of length from I to I + AZ. This VLD obeys a scal- 
ing property n(l) oc l~ a with a — 1.34 ±0.18. Kobayashi 
also found a scaling property of n{l) with a w 1.5 |f|. 

Why can such diverse systems acquire scale invari- 
ant properties? In the case of the large networks such 
as WWW, Barabasi and Albert showed that the scal- 



ing property comes from two generic mechanisms: (i) 
networks expand continuously by the addition of new 
vertices, and (ii) new vertices are more likely to at- 
tach to sites that have more connections than those with 
fewer connections Q. While, there has been no prac- 
tical study that answers this question about the scale 
invariance in quantum turbulence from a generic point 
of view 8] , although many researchers believe that this 
scaling property results from the Richardson cascade pro- 
cess, through which large vortex loops continuously split 
into smaller loops. The aim of this paper is to propose 
an innovative method to study the origin of the scaling 
property in quantum turbulence at zero temperature us- 
ing the Barabasi- Albert model, hereafter the BA-model. 

We now describe the BA-model for large networks. 
Starting with a small number (mo) of vertices, at ev- 
ery time step they added a new vertex with m(< mo) 
edges that link the new vertex to to different vertices al- 
ready present in the system. The probability II that a 
new vertex will be connected to a vertex i was assumed 
to depend linearly on the connectivity fc$ of that vertex, 
so that LT(fcj) = ki/Hjkj. After t time steps, the model 
leads to a random network with t + too vertices and mt 
edges. They calculated the probability P(k) by numeri- 
cal and analytical methods. Both calculations suggested 
that the probability P(k) of this model is proportional 
to fc -3 , which is approximately consistent with the power 
laws observed in various large networks. 

Before applying the above model to quantum turbu- 
lence, we briefly review the nature of quantized vortices 
at zero temperature. First, quantum turbulence can be 
treated as an assembly of vortex loops of various sizes 
when the system is so large that the boundaries have 
negligible influence on the system. Quantized vortices 
are stable, and thus the only ways that the topology of 
vortex loops can change are vortex- vortex reconnections. 
Second, although vortex-vortex reconnections can both 
divide one vortex loop into two (the split type) and com- 
bine two vortex loops into one (the combination type) 
(see Fig. 12 in 0), 

occurrence frequency of the split 
type reconnection is so larger than that of the combina- 



tion type that we can neglect the latter type. Actually, 
this is confirmed numerically for a dilute vortex tangle at 
zero temperature [Io|. Third, the total vortex line length 
is nearly conserved during the dynamics, even when vor- 
tex reconnections occur. The only exceptional event is 
annihilation of small-scale vortices. Reduction of vortex 
line length may also occur due to energy dissipation (e.g., 
phonon emission) caused by the Kelvin wave cascade [llj 
and also by vortex reconnections However, the re- 
ductions in these cases are small relative to the overall 
vortex dynamics that we can neglect the reduction of line 
length 0, 0] . This fact permits us to reach another im- 
portant conclusion: VLD does not change without recon- 
nections. Finally, according to the dynamical scaling law 
of quantized vortices 0] , a vortex loop with a configura- 
tion that is similar to that of another, except that its size 
is reduced by a factor A(< 1), will follow a similar mo- 
tion as the larger loop except the time scale is shortened 
by A 2 . Hence, we can make the following assumption. If 
the scale of the vortex loop 'A' is A times smaller than 
another vortex loop 'B', the probability that the vortex 
loop 'A' reconnects with itself per unit time is (1/A) 2 
times larger than that of 'B\ 

We confirmed the validity of this assumption by di- 
rect numerical simulation of vortex dynamics under the 
vortex filament model as follows. The method of the nu- 
merical simulation is similar to that of previous study 
|l0|. in which we followed the dynamics of vortex fila- 
ments represented by single strings of points. However, 
we adopt here a new algorithm of vortex-vortex recon- 
nection process, based on considerations of crossing lines 
[lil . We prepared a random polygon as a tentative 
initial vortex configuration which follows a probability 
function 

P(ui, • ■ • , Uiv) = Z~ 1 8(u 1 H h u N ) 
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FIG. 1: Schematic of the model. At every time step one, and 
only one, loop is randomly chosen to divide into two daughter 
loops with the probability given by Eq. 10 . If the chosen loop 
has length h, 1ZU and (1 — lZ)k are assigned as the length of 
these two loops. 
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FIG. 2: Solid and slashed lines represent time dependences of 
total line length assuming 11(1) as Eq.10 and eq @, respec- 
tively. 



Here N is the maximum number of points on the initial 
configuration, Uj is the vector from the position of the 
(j — I)-th point on a vortex filament to that of the j- 
th point, S(x) is the <5-function, z is spatial resolution of 
our simulation, s is a scale factor, and the normalization 
factor is given by 

Z = J dux - J dujv(}(ui + Kuat) 
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Because sz corresponds to an average distance between 
neighboring points, we can change the scale of the initial 
vortex loop by controlling s When s > I, this ten- 
tative configuration has larger distances between points 



than the spatial resolution of the simulation. Hence, we 
can complete the initial configuration by interpolating 
points between the tentative points by using an adap- 
tive meshing routine [l(i|. We follow the dynamics from 
the initial configuration and obtain the time t s when the 
vortex makes the first self-reconnection for s = 5 and 
,s = 10. Although the time t s is sensitive to the ini- 
tial configuration, we see that the quantity (tiq) / (t§) 
becomes 3.758 w 4 = 2 2 (i.e., nearly the square of the 
ratio of the size scales), where the brackets represents 
ensemble average of 1,000 samples. This result justifies 
our assumption about the dependence of occurrence fre- 
quency of reconnections on the scale of vortex loops. 

The facts listed above permit us to propose the fol- 
lowing reasonable hypothesis. The observed power law 
of the VLD can be simply understood as a consequence 
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of continuous reconnections, in which the occurrence fre- 
quency is governed statistically by the dynamical scaling 
law of quantized vortices. In this paper, to clarify the 
role of Richardson cascade process, we focus on dilute 
quantum turbulence, in which split type reconnections 
are dominant and correlation of a vortex loop with other 
loops can be negligible. We now show that the decaying 
process of dilute quantum turbulence at zero temperature 
can be described using a simple model similar to the BA- 
model as the following. It must be noted that the system 
we discuss here is different from that concerned by the 
original BA-model. Hence, the two essential points of 
BA-model, increasing the number of the vertices and in- 
equality among the vertices, appear in different ways in 
our model. In the case of decaying process of superfluid 
turbulence, the number of constituents increases by split- 
ting of vortex loops in the system. The inequality among 
the vertices getting linked corresponds to the inequality 
of getting split, which comes from dynamical scaling law, 
in our case of decaying process of superfluid turbulence. 
Starting with a small number (too) of loops, each of which 
have their own length Zj, at every time step we split one 
of them into two daughter loops. When a loop of length 
I is divided into two, the daughter loops have lengths 1ZI 
and (I — 11)1, where 1Z is a random number (0 < 1Z < 1) 
(FigOJ. To incorporate the scaling effect on the occur- 
rence frequency governed by the dynamical scaling law, 
we assume that the probability n that the i-th loop are 
chosen to divide into two depends on its length Zj, so that 
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We call n(Zj) preferential splitting in this paper, because 
II(Zj) corresponds to the preferential attachment intro- 
duced in the original BA-model. We chose the loop to 
be divided using random numbers and II(Z,). When the 
length of a daughter loop was smaller than the cutoff 
length Z m in, we removed the loop from the system. After 
t time steps, the system has an assembly of loops from 
which we can get a VLD n(l). By averaging -/V samp sam- 
ples, we obtained the ensemble average of VLD (n(Z)). 

Unlike the BA-model, which deals with growing net- 
works, it is not easy to analytically derive the power ex- 
ponent a of VLD, due mostly to the super linear prefer- 
ential splitting. Hence, we ran numerical simulations of 
the system. For the simulations, we assumed mo = 1, 
AUnp = 1,000, L = E™ ^ = 10 6 , and Z min = 1. Fig- 
ure |21 shows the resulting time dependence of the total 
line length. The total line length linearly decreases be- 
cause the small-scale loops were removed. Our model did 
not reproduce the observed time dependence of total line 
length |l0j | because the time variable in our model was 
the number of steps in the simulation. Figure |21 shows 
that our model system leads to a power law distribution 
of VLD. The obtained exponent a « 1, which is con- 
sistent with the exponent obtained by Araki et.al. and 
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FIG. 3: Simulated (n{l)) at t = 100,000 (blue open squares) 
and t — 1, 000, 000(red closed circles). The simulations used 
I1(Z) from Eq.@. The slope of the solid line is given by the 
exponent a = 1. 
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FIG. 4: Simulated (n(Z)) at t = 100,000 (blue open squares) 
starting from random < n(l) >(red closed circles) at t = 0. 
The simulations used II(Z) from Eq.©. The slope of the solid 
line is given by the exponent a = 1. 



Kobayashi et.al.. 

Because both the power law and the amplitude of the 
energy spectrum E{k) hardly change during the decay 
H , we expect that a correct model should produce a dis- 
tribution in which the main features are conserved dur- 
ing the decay. Figure |3 clearly demonstrates that VLD 
is independent of time especially in the small Z region, 
indicating that the system self-organizes its stationary 
scale- free VLD. 

Figure 0] shows the obtained VLD that resulted from 
a random initial distribution of VLD. This shows that 
our model does produce a VLD with main features that 
are independent of the initial VLD; that is, this power 
law is very robust. We also confirmed that changing Z m ; n 
hardly affects the power law of VLD as long as the ratio 
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FIG. 5: Simulated (n(l)) at t = 100,000 (blue open squares) 
and t = 1, 000, 000(red closed circles) for the case when 11(7) 
came from Eq.@ (i.e., constant preferential splitting). 

L/lmin is sufficiently large. 

We now address whether or not the model results, par- 
ticularly the stationary power law of VLD, required the 
preferential splitting 11(2) of Eq. ©. To examine this 
question, we investigated how VLD develops without the 
preferential splitting. Instead of Eq. J3J, we assumed 

n ^-4Mr =const - (4) 

By running the same simulations, we found that the re- 
sults lead to a similar power law; however, the exponent 
a, grew monotonically with time (Fig|SJ. 

Our model reproduces the observed exponent well. To 
be accurate, however, our exponent is slightly smaller 
than the measured a « 1.38 ~ 1.5 0, E|. We do not 
understand the origin of this difference, but a likely ex- 
planation is the following. The dynamical scaling law, 
which was introduced through the preferential splitting, 
is a consequence of the localized induction approximation 
of vortex dynamics |l4j . Hence, we implicitly ignored 
the effect of the nonlocal interaction between vortices on 
the dynamical scaling. The nonlocal interaction may be 
needed to get good agreement with measurements. Or, 
we might need to include the combination-type reconnec- 
tion in our model. 

In conclusion, we have proposed an innovative model 
to study the origin of the power law of VLD. The na- 
ture of quantized vortices allowed us to describe the de- 
cay of quantum turbulence with a simple model that was 
similar to the Barabasi- Albert model of large networks. 
The two important generic mechanisms in the Barabasi- 
Albert model, the growth of networks and the prefer- 
ential attachment, corresponded to the continuous self- 
reconnections and the preferential splitting due to the dy- 



namical scaling law, respectively, in our model. Numer- 
ical simulation of our model resulted in a scale-invariant 
distribution of VLD. Furthermore, the exponent was con- 
sistent with measurements that found a VLD approx- 
imately proportional to l^ 1 . Our results suggest that 
only the Richardson cascade and a dynamical scaling 
law arc needed to make the quantum turbulence self- 
organize with a scale- invariant VLD. This means that the 
emergence of power-law VLD is universal phenomenon in 
quantum turbulence beyond Kolmogorov quantum tur- 
bulence. 
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